Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 29 October 2012 (MN style file v2.2) 



Hyper-luminous Reddened Broad-Line Quasars at z ^2 
from the VISTA Hemisphere Survey and WISE All Sky 
Survey 

Manda Banerji^'^*, Richard G. McMahon^'^, Paul C. Hewett^, Eduardo Gonzalez- 



O 

Solares^, Sergey E. Koposov^'"^ 

^ ' ^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CBS OHA, UK. 



^Department of Physics & Astronomy, University College London, Cower Street, London WCIE 6BT, UK. 
^Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge, CB3 OHA, UK. 
' '^Moscow MV Lomonosov State University, Sternberg Astronomical Institute, Moscow 119992, Russia. 



ABSTRACT 



O " 29 October 2012 

^' 

Of 

Q ■ We present the first sample of spcctroscopically confirmed heavily reddened broad- 

^ \ line quasars selected using the new near infra-red VISTA Hemisphere Survey and 

C/3 . WISE All-Sky Survey. Observations of four candidates with [ J — K) > 2.5 and 

K < 16.5 over ~180 deg^, leads to confirmation that two are highly dust-reddened 
broad-line Type 1 quasars at z^2. The typical dust extinctions are Ay ^^2-2. 5 mags. 

■ We measure black- hole masses of ~1O^M0 and extinction corrected bolometric lumi- 
^ ■ nosities of ~10*'' erg/s, making these among the brightest Type 1 quasars currently 

OO . known. Despite this, these quasars lie well below the detection limits of wide-field 

■ optical surveys like the SDSS with iab > 22. We also present WISE photometry at 

3-22^m, for our full sample of spcctroscopically confirmed reddened quasars includ- 

. ing those selected from the UKIDSS Large Area Survey (Banerji et al. 2012a). We 

■ demonstrate that the rest-frame infrared SEDs of these reddened quasars are simi- 
lar to UV-luminous Type 1 quasars with significant hot dust emission and starburst 

■ quasar hosts like Mrk231. The average 12/im flux density of our reddened quasars is 

similar to that of the recently discovered HyLIRG WS'£:i814-f 3412 {z = 2.452) at 
similar redshifts, with two of our reddened quasars also having comparable 22/im flux 
densities to this extreme HyLIRG. These optically faint, heavily reddened broad-line 

' quasars are therefore among the most mid infrared luminous galaxies at z ^ 2, now 

^ . being discovered using WISE. 

Key words: galaxiesiactive, (galaxies:) quasars: emission lines, (galaxies:) quasars: 
general, (galaxies:) quasars: individual 



1 INTRODUCTION 

Dust-reddened broad-line quasars could represent the cru- 
cial missing link between massive ultraluminous starbursts 
identified in far infrared and submillimeter surveys, and UV- 
luminous quasars in optical surveys. The study of these 
rare transitioning systems, has gained new impetus with the 
advent of sensitive wide- field infrar ed surveys such as th e 
UKIDSS Large Area Survey (LAS) (|Lawrence et all 1200 it ). 
VISTA Hemisphere Survey (VHS; McMahon et al. in prep) 
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and Wide Infrared Sur vey Explorer ( WISE) All Sky Sur- 
vey |Wright et al .I2OIOI ). In our previous work (|Banerii et al.l 
(|2012l ): B12 hereafter), we used the UKIDSS LAS to find a 
population of heavily dust-reddened {{J — K) > 2.5 (Vega)) 
quasars at z ~2 corresponding to the main epoch of black- 
hole growth in the Universe. This new sample of 12 con- 
firmed reddened quasars was the first to be assembled using 
a pure near infrared (NIR) selection and pushed 2 mags 
fainter than previous studies using 2MASS, which had se- 
lected q uasars that were considerably less red ((J — K) > 
1.7-2) (|Glikman et al.ll2007l : [c"utri et al.ll200lD 

Sensitive NIR surveys at 1-2/im are well suited to 
searching for broad-line quasars at z ~2, with significant 



2 M. Banerji et al. 



amounts of dust in their host galaxy. At these redshifts, 
the NIR colours sample a portion of the rest-frame spectral 
energy distribution (SED) at ~0.3-0.8/im that is very sen- 
sitive to dust extinction within the quasar host, but not to 
the hot-dust emission f rom the molecular torus at 1-3 /im 
l|Hvland fc AUenl 1 19821 : iNeugebauer et alj Il98/f l. Ground- 
based NIR imaging also allows discrimination of luminous 
Type 1 broad-line quasars that typically appear unresolved, 
and the less luminous Type 2 sources where the host galaxy 
is often resolved. 

Our aim in this work is to use initial data from the new 
VISTA Hemisphere Survey and the WISE All Sky Survey 
completed this year, to devise a colour selection criterion 
for reddened broad-line quasars and subsequently test this 
selection method via spectroscopic follow-up of a small pi- 
lot sample. The infra-red colour selection scheme we use, 
should in principle be effective at isolating dust-reddened 
broad line quasars and the current study offers the first 
practical demonstration of this selection using new data 
from VHS and WISE. In addition, we also investigate the 
rest-frame infra-red SEDs of our full sample of reddened 
((J — K) > 2.5) quasars using data from the WISE All 
Sky Release. Throughout this paper we assume a flat con- 
cordance cosmology with Ho=70 kms~^ Mpc~^, Q,m = 0.3 
and S^A = 0.7. All magnitudes are in the Vega system unless 
otherwise stated. 



2 PHOTOMETRIC DATA 

2.1 VISTA Hemisphere Survey (VHS) 

VHS is a NIR photometric survey being conducted using the 
VISTA telescope in Chile, that aims to survey 18,000 deg^ 
of the southern celestial hemisphere to a depth 30 times 
fainter than the 2MASS survey in at least two wavebands 
(J and K). In the South Galactic Cap, 4500 deg^ will be 
imaged deeper, including the //-band, and will have sup- 
plemental deep grizY imaging provided by the Dark En- 
ergy Survey. The remainder of the high galactic latitude sky 
will be imaged at YJHK and combined with optical pho- 
tometry from the VST- ATLAS survey. In this work, we use 
-180 deg^ data from the VHS VST- ATLAS overlap region 
at 205° <RA< 245° and -9° <DEC< -4.5°. The VHS 
data in this region reaches average 5a /f-band depths of 
18.2 (Vega). 

2.2 WISE 

The WISE satellite has conducted an all-sky survey 
at mid infrared w avelengths between 3.4/xm and 22/im 
(| Wright et all I2OI0I '). WISE specific mid-IR colour selec- 
tion criteria for isolating both populat ions of obscured 
and unobscured AGN have been devised (|Assef et al.ll2Q10l : 
IStern et al1l2012l : lAssef et al.ll2012l '). These authors propose 
a colour-cut of [Wl— W2]>0.8-0.85 with a more conserva- 
tive cut of [Wl— W2]^0.7 suffering from less incompleteness 
but also prone to higher levels of contamination from pop- 
ulations of dusty and normal star-forming galaxies. These 
colour criteria have been shown to be extremely effective in 
selecting AGN at intermediate redshifts but become incom- 
plete when Hof moves into the Wl band at 2; > 3.4. The 



WISE selection identifies AGN independent of their red- 
dening and so needs to be combined with additional colour 
information at shorter wavelengths in order to distinguish 
highly reddened quasars from those that are unobscured at 
UV and optical wavelengths. In addition, morphological in- 
formation from ground-based NIR surveys can be used to 
separate dusty broad-line quasars from Type 2 AGN. 



2.3 This Sample 

Our aim is to test the effectiveness of a joint VHS-f WISE 
colour selection in isolating highly reddened quasars at z~2. 
We therefore identify reddened quasar candidates that sat- 
isfy the following selection criteria: 

• (J-K) > 2.5 and 13< K < 17 (Vega). A bright sample 
with 13 < K ^16.5 (Vega) is prioritized for spectroscopic 
follow-up. The (J — K) colour-cut corresponds roughly to 
E{B -V)> 0.5 at 2 ~ 2. 

• Unresolved or stellar-like point source classifications in 
the VHS NIR images. This corresponds to mergedclass— — l 
in the VHS catalogues. 

• Signal-to- noise of >10 in the /f-band in order to ensure 
reliable morphological classifications. 

• Candidates are required to be red in all the NIR bands 
- i.e. J > H k H > K. 

The above NIR only selection criteria produces ~600 
candidates with K < 17, 447 of which are at K ^ 16.5. Many 
of these sources are expected to be spurious. In order to 
clean the sample, it is matched to the WISE All-Sky source 
catalogue using a matching radius of 3" and requiring a 
>5(T detection in both the WISE 3.4 and 4.6/^m bands. This 
reduces the candidate list to 321 sources. The VHS-f WISE 
detected sources are then required to satisfy the additional 
criterion of [Wl— W2]^0.85, which produces 11 candidate 
reddened quasars with 13 < K < 17, seven of which are 
at 13 < /f ^ 16.5. Finally, the images for all candidates 
are visually inspected to create the sample for spectroscopic 
follow-up. This visual inspection eliminates two sources with 
a close neighbour in the J-band, which may be affecting 
the J-band photometry, leaving five candidates, which are 
summarised in Table [1] 

Before presenting results of our spectroscopic follow-up, 
we first look at the expected effectiveness of our colour se- 
lection criteria presented above. In Figure [T] we show our 
red quasar candidates from VHS+ WISE, with 13 < /S' <17 
in the (J — K) versus (W1-W2) colour-colour plane as well 
as our sample of already confirmed such reddened quasars 
from B12. Also shown are all point-like sources over the same 
area, down to the same flux limit in the /f-band, 142 known 
UV-lu minous quasars with K <17 from the SDSS Data Re- 
lease 7 jSchneider et al.ll2010h . which have been matched to 
VHS and WISE, 114 local ultraluminous infra red galaxies 
from the IRAS IJy sample (|Veilleux et al.ll2002l ') and a sam- 
ple of 73 L a nd T dwarf stars w ith UKIDSS-LAS and WISE 
photometry (|Gelino et al.ll200^ ). We also plot the tracks of a 
UV-luminous quasar, and an evolving elliptical galaxy tem- 
plate modelled as a simple stellar population with zform = 5, 
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Figure 1. (J — K) versus (W1-W2) colour selection of our red quasar candidates. All stellar objects over 180 deg^, are shown in grey. 
The greyscalc represents the density of these objects while the 1% of outliers in the distribution are shown as the individual grey points. 
Known UV-luminous quasars, local Ultraluminous Infrared Galaxies (ULIRGs) from the IRAS IJy sample and known LT dwarf stars 
have also been plotted. We also show the tracks of a typical unreddened quasar as well as an elliptical galaxy template with a formation 
redshift of ^ = 5. Our reddened quasar candidates are marked as the red circles. The small filled circles represent the spectroscopically 
confirmed sample from B12, the two large filled circles are the 2 new confirmed quasars presented in this work and the open circles 
represent all candidates presented in this study down to K <17. 



Table 1. Sample of 13 < < 16.5 and (J - 
quasar candidates from VHS and WISE. 



K) > 2.5 reddened 



Name 




RA 


DEC 


Kvoga 


Redshift 


VHSJ1350- 


-0503 


13:50:37.24 


-05:03:59.4 


15.97 


2.176 


VHSJ1409- 


-0830 


14:09:29.08 


-08:30:58.7 


16.50 


2.300 


VHSJ1504- 


-0442 


15:04:46.25 


-04:42:41.4 


16.50 


No Spec 


VHSJ1518- 


-0501 


15:18:37.85 


-05:01:38.3 


16.13 




VHSJ1543-^ 


-0439 


15:43:21.91 


-04:39:08.6 


16.43 





using an updated version of iBruzual fc Chariot! (j2003l l and 
the publicly available web-based tool EzGAlQ. 

The reddened quasars are significantly redder in terms 
of their (J—K) colours compared to the bulk of point sources 
in the VHS as well as the cool L and T dwarf stars. The red- 
dened quasars are also considerably redder than known UVX 
quasars in the SDSS. Luminous infrared galaxies could en- 
ter our selection box as illustrated by the comparison to the 
IRAS IJy sample in Figure [1] but our morphological classi- 
fication requirement of point sources in the Jf-band means 
that these ULIRGs would have to be extremely compact 
and/or at high redshifts. Finally, compact elliptical galaxies 
too could enter our selection box but only at z > 2 where 
they are unlikely to exist in large numbers and extremely 
unlikely to be brighter than K < 17. 



http:/ /www. baryons.org/czgal/ 



3 SPECTROSCOPIC FOLLOW-UP 

Spectra were taken for four out of our five reddened quasar 
candidates using GNIRS, the NIR spectrograph on Gemini- 
North. We used GNIRS in cross-dispersed mode in order to 
obtain the widest wavelength coverage. Each target was ob- 
served for 20 minutes using 4 exposures of 5 minutes with the 
target dithered along the slit in a classical ABBA pattern 
for purposes of sky subtraction. We used the short camera 
(0.15"/pix) and the 31.7 lines/mm grating with a 1" slit. 
This provided coverage over 0.9-2.5/.tm across six orders. 
Only the two reddest orders are used here corresponding 
to the atmospheric windows at 1.6/xm and 2.2^m. The in- 
strumental resolution is measured to be ~10A or 160 km/s 
from the sky lines. All data were reduced using standard 
GNIRS pipelines (vl.ll) and utilising the Image Reduction 
and Analysis Facility (IRAF) 

Out of the four reddened quasar candidates observed, 
two are confirmed to be broad line quasars at z~2. The spec- 
tra can be seen in Figure [2] along with the best-fit SEDs de- 
rived by fitting to the broadband photometry as in B12. We 
derive dust extinctions, bolometric luminosities and black- 
hole masses for these two quasars using the methods out- 
lined in B12. These are summarised in Tabled These prop- 
erties are fairly typical of our sample of reddened quasars 
selected from the UKIDSS-LAS (B12). Crucially however, 
the VHS targets are in the southern hemisphere and there- 
fore suitable for detailed follow-up using ALMA. Optical 
spectroscopy is needed to confirm the identity of the two 
sources where no lines were seen in the NIR. These could be 
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Figure 2. NIR spectra of our two confirmed dusty broad-line quasars at z~2 from YYIS+WISE. Tlie region of low atmospheric trans- 
mission between the H- and /C-bands is shown as the shaded rectangular area. Vertical lines mark the expected positions of emission 
hues - in order of increasing wavelength - H/3 (486lA), [OIII] 4959A, 5007A and H« (6563A). 



Table 2. Properties of VHS-I- WISE Spectroscopically Confirmed 
Reddened Quasars 





VHSJ1350-0503 


VHSJ1409-0830 


Av (mags) 


2.2 


2.5 


loglo{Lbol/crg s^i) 


47.2 


47.1 


loglo(MBH/M0) 


9.39 


9.22 


Deroddened 


-28.55 


-28.39 


lAB from SED Fit 


22.1 


23.7 



highly compact starbursts with buried AGN, or quasars at 
z < 1.2, 1.75< z <2.00 01 z> 2.8. 



4 MID-INFRARED SEDS OF REDDENED 
BROAD-LINE QUASARS 

We now consider the mid infrared SEDs of all our spectro- 
scopically confirmed reddened quasars including those se- 
lected from the UKIDSS-LAS in B12. In Figure [S] we show 
the SEDs of all 14 spectroscopically confirmed reddened 
quasars in our sample, and comp are them to standard tem- 
plates from PoUetta et al.l (|2007|) including a reddened ver- 
sion of the iPoUetta et al.l (|2007^ QSOl template assuming 
E{B-V) = 0.8 and an SMC-like extinction curve. The SEDs 
are all normalised to a value of 1 at 1 /im, th e well known in- 
flection point in the quasar SED fiElvis ct al.. i, 199411 . The full 
photometric catalogue for all 14 quasars is shown in Table 

El 

We fit two power-laws of the form vFi, oc to the 
rest-frame optical SED at 0.3-lpim as well as the rest-frame 
infrared SED at 1-3/xm. The QSOl template has /3opt = 
0.68 and /3nih, = —0.55. Mrk231 has optical and infrared 
power-law exponents of Popt = —0.91 and /3nir = —0.77, 
while the reddened QSOl template has /3opt = —2.27 and 
/3nir = —1.26. The rest-frame optical power law exponents 
in our sample range from —6 to —0.4 with a median of 
Po-pt ~ —1.47, considerably steeper than the standard QSOl 



template. The different values of Popt reflect the different 
dust extinction values in our sample of reddened quasars. 
The median value of /3nir in our sample is —0.54, consis- 
tent with the standard QSOl template. The infrared slope 
is therefore not particularly sensitive to dust extinction. 
Some quasars have steeper infrared slopes consistent with 
the Mrk231 and reddened QSOl templates. 

The we l l-know n minimum in the quasar SED at ~l/im 
(lElvis et al. 1 1 1994 ) is clearly seen in the rest-frame SEDs, 
and these observations demonstrate the importance of this 
feature on the rest-frame NIR colours. At redshifts of ~0.5 
where the 2 MASS survey has fo und significant numbers of 
red quasars (jGlikman et al.ll2007l ). this minimum falls in the 
J-band therefore making the (J — K) colours of the quasars 
red irrespective of the amount of dust. 

In Figure |4l we also compare the WISE flux densi- 
ties of our reddened quasars to the recently discovered Hy- 
erluminous Infrared G alaxy (HyLIRG), WISE18U+3412 
Eisenhardt et al.ll2012h at a similar redshift (z — 2.452). 
Our quasars are considerably bluer in terms of their 
(W2-W3/W4) colours but the median 12^m flux den- 
sity is 1.46mJy comparable to 1.86mJy in the case of 
W/5'i?1814-|-3412. Two of our most highly reddened quasars 
ULASJ1234-f0907 and ULASJ1539-h0557 (B12) at z = 
2.503 and z — 2.658, also have observed 22^m flux densi- 
ties of >10mJy, comparable to this very luminous HyLIRG. 
These reddened Type 1 quasars are therefore among the 
most mid-infrared luminous sources at z ~ 2, now being 
discovered using the WISE All Sky Survey. 



5 CONCLUSIONS 

We have presented results from a search for reddened broad- 
line quasars using data from the new infrared VISTA Hemi- 
sphere Survey and WISE All Sky Survey. We spectroscopi- 
cally confirm two new reddened quasars in the VHS at z~2 
over an area of 180 deg^. The quasars have broad lines, large 
bolometric luminosities and black-hole masses, and dust ex- 
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Table 3. AB Magnitudes in the UKIDSS YJHK and WISE 3.4, 4.6, 12 and bands for all reddened quasars, along with errors on 

these magnitudes. The full table is published in the online version of this paper. 



Name 


Y 


J 


H 


K 


3.4^m 


4.6^m 


12/im 


22^tm 


ULASJ1234+0907 
ULASJ1539+0557 


26.38±0.84 
21.75±0.36 


22.38±0.29 
21.02±0.20 


20.37±0.21 
19.18±0.08 


18.05±0.03 
17.74±0.02 


17.20±0.03 
16.85±0.03 


16.28±0.04 
16.23±0.03 


14.34±0.04 
14.27±0.03 


13.90±0.17 
13.18±0.06 



Continued online... 
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Figure 3. Normalised SEDs in terms of vF^, of all our reddened 
quasars corrected to the rest-frame and compared to standard 
SED templates from IPoUetta et al.l \200'A \. All SEDs have been 
normalised to a flux density of 1 at 1/im. 
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Figure 4. Rest-frame SEDs of our reddened quasars compared 
to the newly discovered hyp erluminous infra-red gal axy WISE 
1814-1-3412 [z = 2.452) from lEisenhardt et al.l | |2012|) . Our red- 
dened quasars have comparable W3 fluxes to this HyLIRG at 
similar redshifts, with two of our most highly reddened sources 
also showing comparable W4 flux densities. 



tinctions of Av'=2-2.5 mags. We demonstrate that despite 
being intrinsically very luminous, they are too reddened by 
dust to be detected in wide-field optical surveys like SDSS. 
We also present WISE photometry for all our reddened 
broad-line quasars at z~2, including those selected using 
the UKIDSS-LAS (B12), as well as the 0.1-lO^im rest-frame 
SEDs for the entire sample. The rest optical SED power-law 



slopes reflect the different levels of dust extinction in our 
quasars, while the rest infrared power-law slopes are simi- 
lar to Type 1 QSOs with hot dust emission and starburst 
QSO hosts like Mrk231. The reddened quasars are shown 
to have similar 12-22/im flux densities to the recently dis- 
covered hyperluminous infrared galaxy, WISE 1814-1-3412 
at z = 2.452. These reddened quasars are therefore among 
the most mid infrared luminous sources known at these red- 
shifts. Larger samples employing the colour-selection criteria 
presented here, can now be assembled, and complementary 
long-wavelength photometry with Herschel, SCUBA-2 and 
ALMA will shed light on the host galaxy properties of these 
quasars. 
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